On the specific heat of a fermionic atomic cloud in the unitary regime 
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In the unitary regime, when the scattering amplitude greatly exceeds in magnitude the aver- 
age inter-particle separation, and below the critical temperature thermal properties of an atomic 
fermionic cloud are governed by the collective modes, specifically the Bogoliubov- Anderson sound 
modes. The specific heat of an atomic cloud in a elongated trap in particular has a rather com- 
pex temperature dependence, which changes from an exponential behavior at very low tempera- 
tures (T <C Twii), to a T for hu>\\ <C T <C hu± and then continuosly to cx T 4 at temperatures 
just below the critical temperature, when the surface modes play a dominant role. Only the low 
(fiwn >C T <C hu±) and high (hui± <C T < T c ) temperature power laws are well defined. For 
the intermediate temperatures one can introduce at most a gradually increasing with temperature 
exponent. 

PACS numbers: 03.75.Ss 
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Dilute atomic Fermi gases [jj, and especially dilute 
gases interacting with large scattering lengths are under 
intense scrutiny both ex p erimentall y [2, La, |4| and theo- 
retically [E IE 0, H, 0, EE EH IS E3 13, Eli (these lists, 

especially the list of theoretical contributions, is far from 
being exhaustive or, maybe, even representative). While 
quite a number of properties of these systems have been 
clarified, the full and exact picture still awaits to be un- 
covered. What was not appreciated so far is the fact that 
from the information inferred so far, both theoretically 
and experimentally, one can draw a number of unam- 
biguous conclusions concerning the thermal properties of 
these systems. It is widely accepted that the equation 
of state at zero temperature is known theoretically with 
an accuracy of a few percent 0, 0, Q This equation 
of state allows us also to derive the spectrum of low ly- 
ing collective excitations of such systems 0, 0, ^3 1 
namely the sound waves, and one can safely assume that 
experimentally |3| the basic properties of these collec- 
tive modes agree with theory. I shall show here how this 
information can be profitably used in order to extract 
information about the thermal properties of these sys- 
tems at temperatures below the critical temperature for 
transition into a normal state. I shall also show how the 
geometry of the these clouds determines in a somewhat 
unexpected manner the specific heat of these clouds. In 
the first part of this work I shall discuss briefly some of 
the properties of the infinite homogeneous systems and I 
shall turn to the discussion of clouds in harmonic traps 
in the second part. Various aspects of the thermal prop- 
erties of fermionic clouds have been considered by others 
before [l^l EiL Efl]. however, as far as I am aware of, the 
role of the collective excitations below the pairing gap 
has not been discussed previously in this context. 

In an infinite Fermi system with number density 
n = fc^/37r 2 , in the unitary regime, the speed of the 



Bogoliubov- Anderson sound mode is given by |lz 
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where v F = hk F /m and k F are the Fermi velocity and 
wave vector respectively and a is the scattering length. 
The Bogoliubov-Anderson sound waves have the ex- 
pected linear dispersion law with momentum hk, namely 
lu s = ck. The dimensionless parameters ~ 0.44 and 
( s k 1 have been determined recently 0, Ei • At low 
temperatures one can expect that only two types of el- 
ementary excitations exist, the boson-like Bogoliubov- 
Anderson phonons and the fermion-like gapped Bogoli- 
ubov quasi-particles. One can estimate their contribution 
to the total energy E(T) by assuming that at T = the 
system is a Fermi superfluid with a ground state energy 
and a pairing gap determined in Ref. 8]. One thus ob- 
tains [lfij |. see Ref. E3 f° r derivations of such formulas: 
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where e F = h 2 k F /2m, N is the total number of particles 
and the temperature is on the energy scale (ks = 1) and 
A is the pairing gap at T = 0. As illustrated in Fig. I of 
Ref. E3 at T rj Q.2e F the two contributions are essen- 
tially equal (each of the last two terms inside the square 
brackets w 0.06). Off the resonance (when l/k F a ^ 0), 
one has to include corrections similar to those in Eq. QJ. 
Only the (mean-field) exponentially suppressed contribu- 
tion, due to the breaking of the Cooper pairs, and has 
been considered previously 0, 0] . Above the critical 
temperature, such a system behaves as a normal Fermi 
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system and its energy can be estimated with the usual 
textbook formula 



E n (T) = -e F N 
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where w 0.59 can be estimated from the condensa- 
tion energy or be computed directly 0, 0] . Using these 
two expressions, Eqs. and ©, one can show that 
the corresponding free energies cross at a temperature 
Tcross ~ 0.2ef- This fact would suggest that the un- 
known so far critical temperature of a such a system is 
around 0.2ef- This estimate is noticeably lower than 
what a mean-field or BCS-like treatment would suggest 
in the weak coupling limit, namely T c w 0.5ep |5j, [l9j. 
Notice that even in the weak coupling limit the BCS the- 
ory overestimates of the critical temperature |2*oj . and 
the magnitude of this correction seems to be the same 
in the strong coupling limit pj. The above estimate for 
the critical temperature is rather close to a recent ex- 
perimental claim Q . The theoretical mean-field analysis 
performed in conjunction with this experiment [lfij . did 
not consider the role played by the collective modes in 
the thermal properties for T < T c . As it will be shown 
below, in a trap the thermal properties of a fermion cloud 
are modified in a qualitative manner and are expected to 
show the presence of several different regimes. 

In typical experiments with cold fermionic atoms the 
trap is largely harmonic, but it also has a very elongated 
shape. Since tkv\\ <C hui± <C £f, the shape of the atomic 
cloud is that of a very long cigar. The spectrum of the 
collective oscillations of such a superfluid svstems in a 
spherical trap have been evaluated in Refs. 



n, 



W-n(n + Z + 2) + 1, 
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where u> is the frequency of the trap and n and I are 
the radial quantum number and I is the corresponding 
angular momentum. One notices that even though the 
spectrum is distinctly different from a spherical harmonic 
spectrum, nevertheless it is quite close to one (up to the 
notorious scale factor 1 / V3) , namely 
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This behaviour is similar to the change in dispersion from 
lu = vpk for a homogeneous noninteracting Fermi gas to 
u) = VFk/\fi for a homogeneous Fermi superfluid. Since 
I shall be interested in global properties of these systems, 
which depend on the propperties of the entire spectrum 
of these oscillations, the differences between the exact 
and approximate spectra, naively, are not expected to 
change the qualitative picture. For temperatures T <C 
T c = 0(A) = 0(e F ) only such collective oscillations can 
be excited with a noticeable probability. The energy of 



an atomic fermionic elongated cloud is then clearly given 
by the formula 



E s (T)=E(0) + J2 
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tt(n x ,n y ,n z ) = -j= [oj±(n x + n y ) +n z u)\\\ , (7) 

where f3 = l/T and the prime superscript is a reminder 
that only terms with n x + n y + n z > are to be in- 
cluded in the sum. Since the longitudinal wii and radial 
<x>_i_ frequencies differ typically by a factor close to a hun- 
dred, one can distinguish several temperature regimes. 
For very low temperatures, when T <C ^||, the sum in 
the above formula is exponentially small in comparison 
with hu>\\ and the cloud is essentially in the ground state. 
The second regime is for temperatures in the interval 
hjj\ «r« hu>± , when one can easily show that 
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For temperatures in the interval ftuj± <C T <C T c one 
can show that the total energy is essentially quartic in 
temperature, namely 



E S (T) « E(0) + a-. 
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where the constant a can be estimated as the integral 



poo poo poo 

a = 3 3/2 / dx dy dz 
Jo Jo Jo 

3 3 / 2 7T 4 



x + y + z 



30 



= 16.87. 



exp(x + y + z) — 1 
(10) 



Above T c one expects a normal Fermi liquid behavior, 
with an essentially T 2 dependence. Naturally, in this 
rough picture the detailed critical behavior at and around 
the critical point is neglected. The fact that the energy 
as a function of temperature can change its character is 
known in anisotropic solids for a long time |22| , however 
it was not appreciated yet in the context of fermionic 
atomic clouds in anisotropic traps so far. It is important 
to realize that the geometric shape of the could alone can 
lead to this somewhat unexpected behavior of E(T), and 
of the specific heat as well, and that no phase transition 
happens really. The transition from the regime corre- 
sponding to Eq. (JSJ) to the one corresponding to Eq. © 
is not abrupt, but occurs over a rather long temperature 
interval, and can be easily described with a temperature 
dependent exponent. 

The approximate spectrum in Eq. JSJ) does not re- 
produce correctly the dispersion of the surface modes, 
namely f2 sur / rs (An + 3)1/3 (for n <C I) and, as it 
will shown below, their role in the thermal properties 
of a cloud is somewhat unexpected. One can estimate 
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E(T) by evaluating the sum over the quantum numbers 
(n, I, m). For temperatures T 3> one then obtains 



E S (T) = E(0) + J2 



(21 + i)m 



nl 



a. I 



exp(/3M7, i( ) - 1 



£ S M3A0 4 / 3 



142 



2^5 



(ii) 



Both the temperature exponent and the prefactor have 
been obtained numerically. One can evaluate this expres- 
sion also by replacing the sum over I with an integral, 
namely 
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Naturally, there are sub-leading corrections to these esti- 
mates, notably corrections of order 0(T 4 ). It is desirable 
to extend this formula to the case of an axially deformed 
trap. Unfortunately a simple analytical formula for the 
collective modes is not available in this case [2j|. One 
can use the fact that the (local) frequency S} Bur f(S) of 
the surface modes with wave vector k can be determined 
rather accurately from the classical formula 



F(S) 
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where F(S) — |VZ7(r)|g is the force acting on a particle 
at the surface of the cloud, see Refs. for a detailed 
discussion. In the above formula the argument S stands 
for a particular point on the cloud surface. Using an 
axially symmetric harmonic trapping potential U(r) = 
mui 2 (x 2 + y 2 + X 2 z 2 )/2 and the semiclassical formula for 
the contribution of the surface modes to the total energy 
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Here J dS stands for the integral over the surface of the 
cloud. In the case of a spherical trap (A = 1) one recovers 
the previous formula, see Eq. I|12l) . if one includes only 
the modes with n — 0. One can thus conclude that the 
effect of the effect of the deformation of the trapping 
potential can be encapsulated in a simple shape factor. 
Since most experiments are performed in cigar-like traps 
(A <C 1), the role of the surface modes is even larger 
(oc 1/A), when compared to spherical traps. 

It is useful to obtain an estimate of the maximum tem- 
perature up to which this formula is expected to be valid. 



One assumption is that T < A(T), which is likely to hold 
up to temperatures T w T c /2. A reasonable estimate of 
the important angular momenta in the above integral is 
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where l ma x is a rough estimate of the largest single- 
particle angular momentum of a fermion in the ground 
state. Using this estimate one obtains for the total energy 
the following rough estimate 
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Perhaps the thermal energy stored in the surface modes is 
grossly overestimated. At high, but not yet determined, 
temperatures the surface will "melt" and the T 5 regime 
will be replaced by a normal Fermi liquid behavior at 
some point. While a good case could be made that the 
exact numerical factor in front of the second term in the 
above formula is too large, the dependence on the particle 
number is probably correct. One can easily estimate the 
free energy as well and obtain that 

F S {T) = 35^ 



> ^W /3 -35M24iV) 5 / 6 , (17) 

which would have to be contrasted with the free energy 
of the normal phase 
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This estimate of the free energy of the normal phase im- 
plies that no collective modes exist. Since the interaction 
between fermions is attractive, the Landau's zero sound 
mode is unstable. There is still the question of whether 
the first sound mode can exist in such a system at rela- 
tively low temperatures in the normal phase. If the first 
sound modes would be stable, then their contribution to 
the free energy would be very similar to the contribution 
of the Bogoliubov- Anderson sound modes to the free en- 
ergy of the superfluid phase. So far neither theoretical 
nor experimental evidence of a first sound mode in such 
systems in the unitary regime exist. This fact should be 
not construed, however, as evidence of its non-existence. 

If no first sound modes exist, then the inescapable con- 
clusion emerges that the free energy of the superfluid 
phase is likely well below the free energy of the normal 
phase in a trap for a temperature range much larger than 
in the bulk. That would mean that an atomic fermionic 
cloud in the unitary regime remains superfluid at tem- 
peratures larger than the bulk critical temperatures. In 
the above analysis I have obviously neglected the con- 
tribution due to thermally broken Cooper pairs. This 
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contribution, which has been considered previously by 
others, was not the subject of this work. 

In conclusion, the thermal properties of fcrmionic 
atomic gases, especially in the unitary regime, when the 
scattering length greatly exceeds in magnitude the aver- 
age inter-particle separation, are dominated by the ther- 
mal excitation of the Bogoliubov- Anderson sound modes, 
both in the bulk and in traps. In traps, the specific 
heat of such a system, at temperatures below the crit- 
ical transition to a normal phase, has a rather compli- 
cated temperature dependence, ranging from exponen- 
tially damped at very low T's to various power laws. At 
temperatures higher than the axial frequency of an elon- 
gated trap the specific heat has a linear in T regime, due 
to the excitation of predominantly axial sound modes. 
At higher temperatures, close but below the critical tem- 
perature, this regime is replaced with a regime in which 
C oc T 4 , when the excitation of pure surface modes dom- 
inate the thermal properties of atomic clouds in traps. 
This behavior is obviously absent in the bulk and is char- 
acteristic of finite systems only. A similar behavior of 
the specific heat is also expected in the case of an atomic 
Bose gas in a trap, since the spectrum of surface mode is 
similar |2lj . 
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